Coupled-channel pseudo-potential description of the Feshbach resonance in two dimensions 
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We derive pseudo-potentials that describe the scattering between two particles in two spatial dimensions for 
any partial wave m, whose scattering strength is parameterized in terms of the phase shift 8,„. Using our m = 
pseudo-potential, we develop a coupled channel model with 2D zero-range interactions, which describes the 
two-body physics across a Feshbach resonance. Our model predicts the scattering length, the binding energy 
and the "closed channel molecular fraction" of two particles; these observables can be measured in experiments 
on ultracold quasi-2D atomic Bose and Fermi gases with present-day technology. 
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The success in creating Bose-Einstein condensates and 
Fermi degenerate gases has resulted in a renewed interest 
in the scattering between two atoms. At very low tempera- 
tures the de Broglie wavelengths of the atoms are large com- 
pared to the typical Van der Waals length of atom-atom po- 
tentials. As a result, atoms at ultracold temperatures do not 
probe the detailed structure of the interaction potential. Con- 
sequently, the shape-dependent potential can be replaced by 
an appropriate zero-range (ZR) pseudo-potential, which is 
characterized by a single parameter, the generalized scat- 
tering length. For example, s-wave scattering in three spa- 
tial dimensions can in many cases be described accurately 
through Fermi-Huang's pseudo-potential 1 1]. Generalizations 
of Fermi-Huang's pseudo-potential, applicable to higherpar- 
tial wave scattering in 3D, have also been considered (UyllH. 
The foremost advantage of using ZR pseudo-potentials is 
that a number of two-body and some many-body problems 
become analytically solvable, thus highlighting the physical 
meaning of a few key parameters. 

An exciting development in the area of ultracold atom 
physics is the creation of low-dimensional quantum gases. 
The trapping of bosons in quasi-lD geometries, e.g., has al- 
lowed the fermionization of ID Bose gases to be verified ex- 
perimentally |4] and the binding energy of ID molecules to be 
extracted from experiment 1 5]. This paper considers quasi-2D 
systems, for which the motion in the tight confining direction, 
the z-direction, is frozen out llaQl- The equation of state then 
depends on the (generalized) 2D scattering length. A thor- 
ough understanding of the two-particle physics in 2D for any 
partial wave would aid studies of the 2D many-body problem. 
Toward this end, we derive 2D ZR pseudo-potentials for any 
partial wave m and develop a coupled channel model applica- 
ble to two particles with m = across a Feshbach resonance. 

Experiments on ultracold gases routinely utilize Feshbach 
resonances, which allow the interaction strength between two 
particles to be tuned to essentially any value through appli- 
cation of an external magnetic field. A detailed understand- 
ing of Feshbach resonances in 3D underlies studies of, e.g., 
the BEC-BCS crossover and /j-wave pairing. Confinement 
induced resonances in 2D have been predicted 0], and 
recently been observed for p-wave interactions |9], thus al- 
lowing the effective 2D coupling constant to be tuned to es- 
sentially any value including zero and infinity. Assuming a 
strictly 2D geometry, we propose a coupled channel model 



for the lowest partial wave, i.e., m = 0, with four parameters, 
the scattering lengths b\ and b^ of the "open" and "closed" 
channel, the coupling strength p and the detuning e. While 
b\, bi and P are fixed for a given system, varying the de- 
tuning £ corresponds to changing the strength of an external 
magnetic field. To the best of our knowledge, no coupled- 
channel model for 2D resonances has been proposed to date. 
Our model predicts that the dependence of the effective 2D 
scattering length on the external control parameter is distinctly 
different than that of the effective 3D scattering length, thus 
highlighting the non-intuitive behavior of systems with re- 
duced dimensionality. We also determine the binding energy 
and the occupation of the open and closed channel across a 
2D resonance. These observables have been measured across 
a 3D resonance Hot Hill but not yet across a 2D resonance; 
however, we expect measurements on quasi-2D systems to be 
performed in the near future. 

To derive 2D ZR potentials for any partial wave we assume 
that the atom-atom potential V/„ f (p) depends only on the dis- 
tance p between the two atoms, that is, we neglect any angular 
dependence that may arise from spin-dependent interactions. 
The center of mass motion can then be separated off, and the 
radial Schrtidinger equation reads 
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(V(p)-E) 



R m (p) = 0, (1) 



where /j denotes the reduced mass of the two-atom system, 
and V (p) = Vint(p)+Vext(p) [V ex t{p) denotes the p-dependent 
part of an external confining potential, see below; for now we 
set V ex t =0]. In Eq. Q, m denotes the orbital quantum num- 
ber; m = applies to scattering between two 2D bosons or two 
2D fermions with opposite spin, m = 1 to scattering between 
two spin-polarized fermions, and so on. We now derive m- 
dependent 2D ZR pseudo-potentials, V,„;(p) = Vm (p), which 
reproduce the low-energy observables of a shape-dependent 
short-range atom-atom potential. 

We write the pseudo-potential V,,'"(p) in terms of a 8-shell 
of radius s and a yet to be determined operator O m (p) |3], 

V-,f(p) = {<3m(p)5(p-*)W (2) 

The solutions to Eq. (fTJ for V,„ r = Vm can be written in terms 
of the cylindrical Bessel functions J m {kp) and N m (kp). For 
p < s, only J m , which is regular at the origin, contributes, 



R m (p)=B m J m (kp), 



(3) 
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where k — \J2^E /ft 2 . The 8-shell introduces a phase shift 
8 m (£) of the mth partial wave so that the wave function for 
p > s is given by 



R+{p) =A m [J m (kp)-tan(8 m (k))N m (kp)] . 



(4) 



In Eqs. (0 and (|4}, B m and A m denote constants to be deter- 
mined below. 

For s-wave scattering (m — 0), the phase shift So(fc) deter- 
mines the 2D energy-dependent scattering length ao(k) 11211 . 



k 



-cot8 (fc)-y 



(5) 



where y denotes Euler's constant. The energy-dependent scat- 
tering length ao(k), which is always greater or equal to zero, 
is defined such that the scattering wave function has a node 
at p = ao(k). The unusual functional form of ao(k), i.e., the 
exponential dependence on the phase shift, is a direct con- 
sequence of the logarithmic dependence of N m (kp) on kp, 
No(kp) w |[ln(fcp/2) +y] for kp — > 0. For higher partial 
waves, we define generalized energy-dependent "scattering 
lengths" a m (k), which have dimensions of (length) 2 " 1 , as 



.(*) = - 



tan(8,„(Ar)) r(m)r(m + 1)2 
it 2 ™ jT 
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(6) 



Since a\{k) has dimensions of (length) 2 we refer to it as 
scattering area. Energy-independent generalized scattering 
lengths a m are readily defined through a m = lim^o a m (k) . 

Imposing continuity of the wave function R m (p) at p = s, 
that is, requiring R^( s ) = Rm( s )> allows B m to be expressed 
in terms of A m . Integrating the Schrodinger equation from 
p = s — e to s + e, and then taking the limit e — > 0, results in 



ft 2 / 
2/j 



+ / 



(P) 



= O m {s)R m {s). (7) 



Plugging Eqs. (0 and @ into Eq. Q and taking ks <C 1, de- 
termines the operator O m (s), and hence the pseudo-potential 
Vm J (p). For m = 0, we find 



v ps (p) 



-£tan(5o(*)) 



l-^^Mk,p))np d P 



p8(p - s) 



,(8) 



where fo(k,p) = 1 +y+ln(A;p/2). The explicit k dependence 
of the m = ZR potential drops out when Vg is written in 
terms ao(k). Our s-wave pseudo-potential agrees with the A 
potential derived in Ref. 11311 if one sets A equal to k (see also 
Ref. I14ln . Below, we use the boundary condition implied 
by the pseudo-potential Vq" to develop a coupled-channel ZR 
model for 2D scattering. 

For m > 0, a straightforward yet somewhat tedious calcula- 
tion gives 



ft 2 [r(m+l)] 2 tan(8 m (A:)) (2\ 2m 



(in a 2 - 
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tan(s " ,w) / m (fc,p))p m+ ' 3 p 2 "' p 5(p 5) 
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where 



f m (k,p) =\(/(m)-21n 



kp 



r=2m— 1 



2m — r 



(10) 



Here, vff(m) = V|/(l) +\\t(m+ 1), where \\t denotes the digamma 
function. As written in Eq. the m > pseudo-potential 
leads, despite the regularization operator, to divergencies at 
p — > if the energy-dependent coefficients C m and D m of the 
p~ m and ln(p)p" ! terms in the expansion of the eigenfunction 
sought differ from the corresponding coefficients c m and d m 
of the expansion of the irregular free particle solution N m (kp). 
To cure this divergence, the right hand side of Eq. (|9jl has to be 
multiplied by d m C m /(c m D,„), resulting in a pseudo-potential 
that has to be determined self-consistently 1 15] . A similar 
self-consistency condition is not needed for systems with odd 
dimensionality. 

The bound state energies E^" d of the pseudo-potentials Vm 
can be determined through analytic continuation 1 3|. For m = 
0, we recover the well known expression for the ZR binding 

H, 



energy E l ^ lnd 



-'bind 



-h 2 



4exp(-2y). 



2n(a (k)) 2 

For m > 0, the binding energies E l ™ ld are given by 
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r(m)r(m + i) 



1/m 



(11) 



(12) 
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The occurance of an ";" in Eq. \\2\ appears odd at first sight. 
This puzzle is resolved by noting that the 2D generalized 
energy-dependent scattering lengths a m (k) with m > are, 
in contrast to the ID and 3D counterparts, complex. Conse- 
quently, the binding energies Ef" ld are determined by those 
a m (k) for which the real part vanishes. For m = 1, e.g., 
we find excellent agreement between the exact binding en- 
ergy determined from the energy-dependent scattering area 
via Eq. (I12> and that for a square well potential with range 
Po> Vi n t(p) — Vsw{p) = —Vo for p < Po and zero otherwise. 

We now use the proposed pseudo-potentials to determine 
the eigenspectrum of two atoms in 2D under external har- 
monic confinement, that is, in Eq. Q we consider V ext (p) = 
/j<D 2 p 2 /2. Such a two-body system can be realized experi- 
mentally with the aid of a ID optical lattice with doubly- 
occupied lattice sites, for which the tunneling between neigh- 
boring sites is neglegible. The solutions for p < s and p > s 
are proportional to the confluent hypergeometric functions M 
and U, i.e., R m (p) - p» exp(-p 2 /2a 2 )M(-x, m + 1 , p 2 / fl 2 J 
andi?+( p)oc p m exp(-p 2 /2a 2 )t/(-x,OT+l,p 2 /a 2 J, where 
flfto = y ft/ (/jco) denotes the oscillator length and % a non- 
integer quantum number, E nfjL = (2% + 1 +m)ftco. Enforcing 
continuity at p = s determines the eigenenergies E mx implic- 
itly in terms of a m (k), 
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not, as done in 3D to be proportional to the wave func- 
tion components themselves. The scattering length a£ (E,e) 
predicted by Eqs. fl!4i and d!5t in the limit co — > is 



(£,e) = foiexp -p 2 / 



Y--in( fjf(e-£) 



•(16) 



To determine the behavior of the scattering length Oq c (E , e) 
in the vicinity of the resonance as a function of the mag- 
netic field strength B, we Taylor-expand Oq C (E,e) about the 
resonance position Er, which is given by the binding en- 
ergy E^ nd of the strongly closed molecular channel, Er = 
2exp(— 2y)h 2 / (fjb^)- We find a simple functional form for 
Figure 1: E l% as a function of a { (k = 0) (co = 2jw). Asterisks show a% c (E = 0,5) in terms of the background scattering length 

Eq- 1131 . for the energy-dependent pseudo-potential V/ M (p) and 
solid lines show E\-^ for V$w with po = 0.01a/, o . For comparison, 
horizontal dotted lines indicate the energy levels for a\ = 0. 
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Aj,„, the resonance width A and the resonance position Br, 



(E = 0,B) =A bg exp[-A/(B-B R )], 



(17) 



where k,„^ = ^/2/jE,„ x /h and (x) r = x(x + 1) — (jc + r — 1) 
with (x)o = 1. Equation (1131 contains the generalized energy- 
dependent 1 17] 2D scattering length a m (k) evaluated at k = 
k nfl and is valid for m > 0; for m — the eigenenergies are 
given by Eq. dl8t with (3 = and b\ — ao(k). 

Asterisks in Fig. \l\ show the m — 1 eigenenergies for 
two particles interacting through the proposed ZR potential, 
Eq. Jl 31 . under harmonic confinement over a large range of 
zero-energy scattering areas a\. For comparison, solid lines 
show the exact eigenenergies, determined semi-analytically, 
for two particles under harmonic confinement interacting 
through a square well potential with range po = 0.01a/, o . Fig- 
ure [2 illustrates excellent agreement between the eigenen- 
ergies for two 2D particles interacting through the pseudo- 
potential and those interacting through the square well poten- 
tial. We find similar behaviors for higher partial waves. 

To describe the two-body physics across a 2D Feshbach res- 
onance for two bosons or for two fermions with opposite spin, 
we develop a coupled-channel ZR model for m = 0. For p > 0, 

the wave function with components R^\p) and (p) satis- 
fies the harmonic oscillator Schrodinger equation in the rela- 
tive coordinate, 



if # 

u \ dp 1 



1 9 9 



file 3 2 



I 

P dp 



I 3 

^-iuco 2 p 2 



R 



= 0, (14) 



where e denotes a detuning (e > 0), which can be changed, 
e.g., by varying the strength of an external magnetic field. The 
2D pseudo-potential Vq S imposes a boundary condition at p = 
0, which we parameterize as 



HP/bi) P 



Ps^o 

r, 3 p(2) 




(15) 



The coupling between the two channels is characterized by 
the dimensionless parameter p. To ensure a divergence-free 
treatment, Eq. ( 1151 takes this coupling parameter to be propor- 
tional to the derivative of the wave function components, and 



where Ai, g =b\, B R = Er/M and A = 2$ 2 Er/M (M denotes 
the difference in magnetic moment between atoms in the open 
and closed channel). Just as the 2D single-channel scatter- 
ing length ciq [Eq. Q], the functional dependence of the 2D 
coupled-channel scattering length Oq C on B near a resonance 
is distinctly different from that of the 3D counterpart. In prin- 
ciple, the parameters b\, bi and p can be determined by com- 
paring Eq. dl7> with experimental data for a specific 2D Fes- 
hbach resonance. Since no such data exist to date, the inset 
of Fig. |2 illustrates the behavior of a£ c as a function of £ for 
E = 0, a\ = 0.5o/, o , a? = 0.05fl/,„ and p = 0.1. The scatter- 
ing length Oq C changes from infinity to zero at the resonance 
value (indicated by a vertical dotted line). 

We find an implicit eigenequation for the eigenenergies Eq„, 
Eq„ < £, of the coupled-channel ZR model for two particles 
under harmonic confinement with m = [Eqs. dl4> and J15H . 
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Figure [5] shows the eigenenergies Eq„ as a function of £ for 
b\ = 0.5a/ 7O , £>2 = 0.05a/,„ and three different values of p, i.e., 
P = 0.03 (solid lines), 0.1 (dotted lines) and P = 0.3 (dashed 
lines). The non-vanishing coupling leads to a series of en- 
ergy level crossings at £ w Er, which become narrower as 
the coupling p decreases. For higher-lying states, i.e., larger 
Eq„, the resonance of the energy-dependent scattering length 
% (e, E), Eq. d!6i . moves to larger £ values, and thus explains 
why the energy level crossings move to larger £ for higher n. 
The energy of the state n — 1 for E < Er corresponds to the 
binding energy of the open channel, i.e., of channel (1). We 
checked that a 2D coupled channel square-well system (see 
Ref. l2(ill for the 3D analog) reproduces the results obtained 
for our proposed coupled-channel ZR system. 

Our coupled channel model leads to a mixing of the 

(2) 

strongly closed molecular channel R J and the open channel 

RqJ . To quantify the admixture of the strongly closed molecu- 
lar level across the resonance, we define the molecular fraction 
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Figure 2: Molecular fractions />o (solid line) and Pi (dotted lines), 
Eq. 1191 . for b\ = 0.5a/, o , £>2 = 0.05a/, o and different coupling con- 
stants p as a function of e. Inset: 2D scattering length Og C (£,e) as a 
function of e for E = 0, b\ = 0.5a/, o , hi = 0.05a/, o and (3 = 0.1. 



The main part of Fig.[2]shows P„ for b\ = 0.5a/j O , bj_ = 0.05a/, o 
and different p values as a function of e; a solid line shows Po 
for p = 0.1, and dotted lines show Pi for p = 0.03, 0.1 and 
0.3. The molecular fraction Po is close to one for small e, and 
drops to zero as e is swept across the resonance. This indicates 
that the character of the n = state changes from "strongly 
closed molecular" to "weakly closed molecular" as e changes 
from e < e« to e > Zr. The molecular fraction of the n = 1 
state is close to zero away from resonance for all coupling 
strengths considered. Near resonance, however, Pi depends 
on (3. For weak coupling, Pi approaches one on resonance. 
For strong coupling, in contrast, Pi is comparatively small as 
shown in Fig.[2]for p = 0.3. This suggests that the 2D ana- 
log of the BEC-BCS crossover can be best studied utilizing 
2D Feshbach resonances with strong coupling for which the 
admixture of the strongly closed molecular channel is small. 
Similar behavior is found for 3D systems fl9ll . 
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Figure 3: Eigenenergies Eq„, Eq. 1 1 8i . for b\ = 0.5a/, o , bi = 0.05«/ ]o 
and three different values of P, i.e., P = 0.03 (solid lines), 0. 1 (dotted 
lines) and 0.3 (dashed lines) as a function of £. 



In summary, this paper derives a series of 2D pseudo- 
potentials, which describe the low-energy scattering of two 
particles with partial wave m. The boundary condition im- 
plied by the m = pseudo-potential is then used to develop 
an analytically solvable coupled-channel model, which de- 
scribes the physics across a 2D Feshbach resonance. The pre- 
dicted dependence of the effective 2D scattering length Aq C 
on the strength B of the external magnetic field, Eq. il/l . 
may prove useful in analyzing experimental data for quasi- 
2D Bose gases or two-component Fermi gases. We also deter- 
mine the binding energy of 2D dimers, which can be measured 
with present-day technology by utilizing optical lattices with 
doubly-occupied lattice sites in the no-tunneling regime 1 1 Oil . 
across a Feshbach resonance. 
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